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Abstract 


Neuronal modeling of patch-clamp data is based on approximations which are valid under specific assumptions regarding cell 
properties and morphology. Certain cells, which show a biexponential capacitance transient decay, can be modeled with a two- 
compartment model. However, for parameter-extraction in such a model, approximations are required regarding the relative sizes of 
the various model parameters. These approximations apply to certain cell types or experimental conditions and are not valid in the 
general case. In this paper, we present a general method for the extraction of the parameters in a two-compartment model without 
assumptions regarding the relative size of the parameters. All the passive electrical parameters of the two-compartment model are 
derived in terms of the available experimental data. The experimental data is obtained from a DC measurement (where the 
command potential is a hyperpolarizing DC voltage) and an AC measurement (where the command potential is a sinusoidal 
stimulus on a hyperpolarized DC potential) performed on the cell under test. Computer simulations are performed with a circuit 
simulator, XSPICE, to observe the effects of varying the two-compartment model parameters on the capacitive transients of the 
current response. Our general solution for the parameter-estimation of a two-compartment model may be used to model any 
neuron, which has a biexponential capacitive current decay. In addition, our model avoids the need for simplifying and perhaps 
erroneous approximations. Our equations may be easily implemented in hardware/software compensation schemes to correct the 
recorded currents for any series resistance or capacitive transient errors. Our general solution reduces to the results of previous 
researchers under their approximations. © 2002 Elsevier Science B.V. All rights reserved. 
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Nomenclature 

T single-compartment time constant (s) 

TO CpRcRp/(Rc + Rp) (s) 

ra [((1/Cp)(/Ruy + 1/Rco+1/Rs)]~ | (s) 

T see Eq. (14) (s) 

T3 two-compartment time constant (slow) (s) 
T4 two-compartment time constant (fast) (s) 
Rs pipette resistance (ohm) 

Rui resistance of compartment M (ohm) 

Cu capacitance of compartment M (farad) 
Rp resistance of compartment D (ohm) 

Cp capacitance of compartment D (farad) 
Re resistance connecting compartments M and D (ohm) 
Vo command potential (volts) 
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Vi(t) potential at compartment M (volts) 
V3(t) potential at compartment D (volts) 
Vp steady-state voltage component of V(t) (s) 
ii(t) total input current through the pipette (ampere) 
iss steady-state component of i,(t) (ampere) 
Ao single-compartment exponential-decay coefficient (ampere) 
A, two-compartment current exponential-decay (slow) coefficient (ampere) 
A> two-compartment current exponential-decay (fast) coefficient (ampere) 
é see Eq. (28) (s~”) 
7) operating frequency of the sinusoidal stimulus (radians/s) 
Y(@) single-compartment admittance at a frequency w (ohm ') 
Ar(@) real part of Y(w) (ohm 7 ') 
Br(@) imaginary part of Y(@) (ohm ~') 
RrorTaL single-compartment thevenin equivalent resistance (ohm) 
Rrotr two-compartment thevenin equivalent resistance (ohm) 
Yror(@) two-compartment admittance at a frequency w (ohm~') 
Zror(@) two-compartment impedance at a frequency m (ohm) 
a(@) real part of Yror(@) (ohm~') 
B(@) imaginary part of Yror(@) (ohm ') 


1. Introduction 


In most electrophysiological studies, whole-cell vol- 
tage-clamp techniques are an established method to 
measure ionic currents from single cells. An ideal 
voltage clamp has two functions: first, it imposes a 
command potential on the membrane such that the cell 
membrane potential is equal to the command potential 
and second, it measures the ionic current. In reality, 
there is considerable resistance to the flow of current 
through the pipette to the cytoplasm. As a result of this 
current flow, the membrane potential is not equal to the 
applied command potential. The signal response time 
(time constant) is essentially a function of the pipette 
resistance, the membrane resistance and capacitance. An 
estimation of these electrical components is needed to 
compensate for the series resistance and capacitive 
transient errors in a voltage-clamp experiment. 

Compartmental models are often used to represent a 
cell’s passive electrical properties, which provide a useful 
insight into the different cell conduction mechanisms. In 
the simplest case, a single-compartment model (Fig. 1) 
can be used to represent the electrical characteristics of 
small cells (longest dimension <100 um). Here, the 
capacitive transient of the current response to a hyper- 
polarized voltage step (i.e. a negative step voltage 
relative to the intracellular medium) has a single- 
exponential decay time constant. Also, the membrane 
resistance is usually 50-100 times larger than the pipette 
resistance (5—10 MQ). As such, the method of para- 
meter-estimation in a single-compartment model be- 
comes quite trivial and, over time, numerous methods 
have evolved for the extraction of these parameters. If 
we assume the membrane resistance is much larger than 


the pipette resistance, then all single-compartment 
model parameters can be estimated within reasonable 
accuracy (Sakmann and Neher, 1994). A more exact 
method of parameter-estimation involves admittance 
measurements (Lindau and Neher, 1988), where the 
model parameters can be expressed in terms of the real 
and imaginary components of the measured admittance. 
Also, the dynamic changes in passive cellular character- 
istics may be measured with phase-tracking techniques 
(Fidler and Fernandez, 1989) or by other dual-frequency 
measurement approaches (Donnelly, 1994). Computer 
simulations have been performed by Sala et al. to 
examine the sources of errors introduced by the single- 
compartment model parameters in single-electrode vol- 
tage-clamp experiments (Sala and Sala, 1994). Recently 
a method, which employs planar pipettes in silicon, was 
suggested to reduce noise sources in patch-clamping and 
for high-throughput screening (Pandey and White, 
2001). 

A single-compartment model cannot describe some 
cells. In neurons, such as the cerebellar Purkinje cells, 
the capacitive current response to a hyperpolarizing 
voltage pulse is a sum of two exponentials. These cells 
can be modeled satisfactorily with a two-compartment 
equivalent circuit. In a two-compartment model, the 
equivalent electrical circuit consists of two main com- 
partments (soma and proximal dendrites on the one 
hand, distal dendrites on the other) connected by a 
resistance. In hippocampal pyramidal cells, three ex- 
ponentials are needed to model the capacitive decay. In 
these instances, and many more complex cases, a cable 
analysis may be more appropriate than a multi-com- 
partment model (Koch and Segev, 1989). The cable 
method for analyzing passive electrical data from 
neurons consists of decomposing the voltage response 
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Fig. 1. (A) A single-compartment model showing the passive electrical 
equivalent circuit of a small cell. (B) The typical waveforms of the 
command voltage Vo and the recorded current i,(t) versus time. 


of the cell into a series of exponential functions and 
substituting the time constants of these exponential 
functions into equations derived from cable theory 
(Rall et al., 1992). As an alternative to this traditional 
method of cable analysis, the method of ‘integrals of 
transients’ has been employed to analyze the passive 
electrical data from neurons (Engelhardt et al., 1998). 
Yet, the simpler the compartmental model, the easier 
it becomes to estimate its model parameters (Roth and 
Hausser, 2001). Unfortunately, even in a two-compart- 
ment model, some approximations need to be made to 
simplify the model and thereby, ease the method of 
parameter-extraction. Llano et al. (1991) employed a 
two-compartment model to fit the capacitive transients 
of the excitatory currents from Purkinje cells as a sum of 
two exponentials. Llano et al. (1991) assumed the 
resistances in the proximal and distal compartments 
were very large compared with the other resistances in 
the equivalent circuit. The biexponential decay of 
capacitive transients from hippocampal cells was de- 
scribed using a two-compartment model by Mennerick 
et al. (1995) assuming the resistance of the distal 
compartment is very large. In another work, a simplified 
two-compartment model was used to extract the mem- 
brane properties of bipolar neurons, assuming very high 
resistances in the proximal and distal compartments 
(Mennerick et al., 1997). A recent work compared the 
various approximations made by previous authors while 
using a two-compartment model, summarizing the 


validity and limits of these approximations under 
different conditions (Nadeau and Lester, 2000). 

The reliability of any parameter-extraction technique 
depends on the approximations made while extracting 
the model parameters and on the validation of such 
approximations for the given cell and the bandwidth of 
operation. Compensation circuitry, either included in 
the software/hardware of the patch-clamp setup, relies 
on the accuracy of these parameter-extraction techni- 
ques to compensate for the errors introduced by the 
circuit resistance and capacitance (Traynelis, 1998). In 
some cases, certain simplifying approximations regard- 
ing the size of the model parameters can be justified. But 
in other cases, especially when dealing with complex and 
larger cells, a more exact model for the estimation of the 
parameters is required. This points out the need for a 
general method of parameter-estimation in a_ two- 
compartment model, which is not subjected to over- 
simplifying approximations. 

In this work, we present an exact, analytical method 
for the parameter-estimation of a two-compartment 
model. All the passive electrical parameters of the two- 
compartment model are derived based on the DC 
measurements (where the command potential is a DC 
hyperpolarizing voltage) and on AC measurements 
(where the command potential is a sinusoidal stimulus 
resting upon a DC hyperpolarizing voltage). Our 
equations are compared to those derived by previous 
authors and, under approximations used by these 
authors, we obtain the same results. Computer simula- 
tions are performed on a circuit simulator, XSPICE, to 
demonstrate the contributions of varying each model 
parameter on the capacitive transients of the current 
response. 


2. Single-compartment model (a review) 


We assume a reasonably small cell (< 100 um) with a 
membrane resistance Ry, and a membrane capacitance 
Cm. The command voltage Vo is applied to the cell 
through a pipette with a series resistance of Rs. Fig. | 
shows the electrical equivalent circuit of the small cell, 
along with the voltage and current waveforms asso- 
ciated with the circuit. This equivalent circuit represen- 
tation of a small cell is often used with patch-clamp 
techniques in electrophysiological experiments. We 
assume the cell membrane is isopotential and there are 
no voltage-dependent, active conductances. In small 
cells, the typical values of the single-compartment model 
parameters are Rs = 10 MQ, Ry = 1 GQ and Cy = 15- 
100 pF. In such cases, the assumption of Ry > > Rg is 
valid and simplifies the model. We will describe briefly 
the standard equations of this single-compartment 
model, estimating the model parameters by making the 
assumption of Ry > > Rs. Next, we will review a more 
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accurate method for the estimation of the single- 
compartment model parameters based on admittance 
measurements with no simplifying assumptions. We will 
extend this approach to a two-compartment model, 
which will provide a technique for the extraction of the 
two compartment model parameters. From Fig. 1(A), 
we derive the equations 


vO=5% (1 exp ")) (1) 
S\*M 


is the circuit time constant and the current is 


Vi t 
i,(1) = 2 (1 ! ) +107 exp (-*) 2) 
Rg RsCy)  R2Cuy t 
The experimentally recorded current i,(t) has the form 
i,(t) = igg + Ag exp(—t/t) (3) 


where iss is the steady state component of ij(f). 
Comparing Eq. (1) and Eq. (3) we have 


V, 
en ee and 4, =—& (4) 


With the assumption Ry > > Rs, the model parameters 
can be written as (Sakmann and Neher, 1994) 
Ry=—2, Rez 2, Cy as (5) 
Iss i,(0) Rs 

Computer simulations of the single-compartment 
equivalent circuit are performed on XSPICE, a circuit 
simulator. Fig. 2(A) shows the effect of varying the 
series resistance on the capacitive transient of the 
current i,(¢) under the application of a 10 mV command 
voltage step. The values of the model parameters are 
chosen as Ry = 100 MQ, Cy =20 pF and Rg =2, 3.5 
and 5 MQ. We see the series resistance compensation is 
increased (i.e. as Rs is decreased), the time constant t 
decreases and the value of i,(0) or Ao increases 
(according to Eq. (2)). In Fig. 2(B), the amplitude of 
the command voltage step is varied and its effect on the 
capacitive transient of the current i;(¢) is observed. The 
values of the model parameters are Ry =100 MQ, 
Cm =20 pF and Rg =2 MQ. In Fig. 2(B), we observe 
that even though the time constant is fixed (t = 39 us), 
the current i,(0) or Ao increases with increasing 
magnitude of the command voltage Vo (according to 
Eq. (2)). 

We now review a method of parameter estimation for 
a single-compartment model, which gives an exact 
solution for the model parameters, without simplifying 
assumptions (Lindau and Neher, 1988). In this techni- 


nA 


(B) 


Fig. 2. (A) The effect of series resistance compensation on the 
capacitive transient of the current i,(¢). The model parameters are 
Ry = 100 MQ, Cy = 20 pF and Rs = 2, 3.5 and 5 MQ. The fitting 
parameters for the single-exponential decay of the capacitive transients 
are: (a) for Rs = 2 MQ, t = 39 ps, Ag = 4.9 nA, igg = 0.098 nA; (b) for 
Rs = 3.5 MQ, t= 68 ps, Ao = 2.7 nA, iss = 0.096 nA; (c) for Rs = 5 
MQ, t= 96 ps, An = 1.9 nA, igg = 0.095 nA. (B) The effect of varying 
the command step voltage Vo on the capacitive transient of the current 
i,(t). The model parameters are Ry = 100 MQ, Cy = 20 pF and Rs = 
2 MQ. The fitting parameters for the single-exponential decay of the 
capacitive transients are: t = 39 pts (a) for Vo = —10 mV, Ap = 4.9 nA, 
iss = 0.098 nA; (b) for Vo = —30 mV, Ao = 14.7 nA, iss = 0.29 nA; (c) 
for Vo = —60 mV, Ap = 29.4 nA, isg = 0.59 nA. 


que, a sinusoidal voltage stimulus resting upon a 
hyperpolarized DC potential is applied as the command 
voltage. The magnitude and phase-shift of the resulting 
current sinusoid are analyzed with a phase-sensitive 
detector to give the real and imaginary current compo- 
nents. These current components divided by the stimu- 
lus voltage amplitude give the real and imaginary 
admittance values. Using this information, along with 
the information from the measured DC current, all the 
parameters of a single-compartment model can be 
exactly determined. 
From Fig. 1(A), the net admittance is 


Y(@) = Ap(@) + jBR(@) 
_ (1 + @ Ry RpCu) + JORM CM 
RU + w@? REC) 


(6) 
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where 
Ree eRe and fo (7) 
Ipc Ry + Rs 
With the experimental values of Ry, Ap and Br, Eq. (6) 
and Eq. (7) give the model parameters as 
Ag — Ry! 
AR + BR — AgRg!’ 
Ry{(dp — Rz'Y + BR} 
AR + BR — AgR;! 
Ge ( I ) (AR + Ba - Ag Rr! 
Oc Br] (Ag — Ry'Y + BR 


Rs 


Ru 


? 


(8) 


where fco=@c/2n is the frequency of the applied 
sinusoid. 

To estimate the three parameters Rs, Ry and Cy of 
the single-compartment model, at least three equations 
are needed. In the Lindau—Neher technique, the admit- 
tance measurement provides two equations and the DC 
current measurement provides the third equation needed 
to estimate all the model parameters. We use a similar 
methodology in the parameter-extraction of the two- 
compartment model. The current response i;(t) to a DC 
command voltage, together with the admittance mea- 
surements, will provide all the information needed to 
exactly estimate all the six parameters of the two- 
compartment model, as shown in the next section. 


3. Two compartment model (DC measurements) 


In this analysis, we assume a cell that can be modeled 
by an electrical equivalent circuit shown in Fig. 3. The 
DC step command voltage Vo is applied to the cell 
through a pipette with a series resistance Rg. The actual 
cell can be represented by two compartments: compart- 
ment M (for the soma and proximal dendrites) and 
compartment D (for the distal dendrites). A series 
resistance Rc connects the two compartments. Assum- 


Fig. 3. Two-compartment model showing the two compartments: 
compartment M (for the soma and proximal dendrites) and compart- 
ment D (for the distal dendrites) connected by a resistance Rc. The 
pipette resistance is Rs and the command voltage stimulus is Vo. 


ing the pipette capacitance is fully compensated, the 
following equations can be derived for the two-compart- 
ment model: 

iy =(Vo- V)/Rs, i, =(V, —V)/Re, l= V,/Ru, 
is =V,/Rp, ig =CydV,/dt, ig = CpdV,/dt, 
W-h=kt+h, h=is tie 


(9) 


From the above set of equations, we obtain two first- 
order differential equations 


dV (t t . ta Vi) V, 

Ce + v0 ee )- f(D ° (10) 
dt Ry Re Re Ra” Te 
dV,(t t. a Vit 

. + vs0( + )- 10) (11) 
dt Re Rp Ro 


Eliminating V,(t) from these two differential equations 
gives a second-order non-homogeneous differential 
equation in V3(t) as 


2 
PHO MOP yof 1 | 


dr? di lt. 3 tty Cy uke 
V, 

= 2 (12) 
Cu Cp Rc Rs 

where 

tJ ( ) 

T Cy \Re Rp 
1 1 1 1 1 

and — ( +—-+ ) (13) 
tT Cy \Ru Re Rs 


The solution of such a differential equation is of the 
form: V(t) =c, exp(—t/t3) +c exp(—t/t4)+ Vp, where 
Vp is the particular solution of the above differential 
equation and the exponential terms are the solutions of 
the homogeneous differential equation. The constants c 
and c, can be determined from the given initial condi- 
tions. For Eq. (12), the particular solution Vp is 


_ VoRmut 
Co(Ry Rs + Ry Re + Re Rs) 


Vp 


where 


=1 
(Rs + Rw) aa 


T = C| + 
Ry Ry Rs + Ru Re + RcRs 


The homogeneous solution of the differential Eq. (12) 
gives the two circuit time constants as 


1 1 1 1 1 

T3 Tot; RECyCp To 7 

( ) ( M/ 
— + - 

t, To Ty Tot, RECyCp 


x (+) (16) 
To TY 
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We also have a relation between the two time constants 
as 


1 1 1 1 
Gaon wm 
Ty «3 To T 
The initial conditions of V2(0) = 0 and V;(0) = 0 give the 
general solutions of V,(t) and V3(t) as 


V(t) = ( VoRutr ) 
; Cp(RyRc + Ry Rs + RRs) 


a ce 
Cana) Hp 


V(t) = ( VoRuRety ) 
To(Ru Re + Ru Rs + RRs) 


x ji+( “2)exp( -) 
Ty — 13 T3 
(= _ “oe ( ‘) (19) 
Ty — T3 T4 


In the approximation of Cy =0 (Nadeau and Lester, 
2000), we observe the time constants in Eq. (15) and Eq. 
(16) reduce to 13 = 12 and t4=0. If we place the values 
of these time constants into Eq. (18), then we will have 
the same result obtained by Nadeau and Lester (2000). 
The next section will provide additional validation of 
the above set of equations for the two-compartment 
model. Now, the experimentally recorded current i;(t) 
has the form 


i,(t) = igg + A, EXP (-£) + A, exp (-£) (20) 


T3 T4 


where iss is the steady-state current component of i,(f), 
while 73 (slow) and ty (fast) are the two time constants of 
the biexponential capacitive transient. Fitting the ex- 
perimental plot of i,(t) with Eq. (20) will provide the 
values of all the equation unknowns (igs, 73, 74, 41, 42). 
From the experimental plot of i,(¢), the following 
parameters can be determined: 


ist Vo R= Vo _ Vo 
oo eg : 1,0) Ay + Ay + fgg’ (21) 
Vo 1,74(A, + Ay + igs)” 


Cu 2/3 = 
R(di,(0)/dt) Vo(A,t4 + A573) 


Eq. (21) allows us to estimate the parameters Rs and Cy 
from the known values of the current parameters (igs, 73, 
t4, A}, Ar) in Eq. (20). Also, from Fig. 3 and Eq. (21), 
Ry can be expressed in terms of Rrot, Rs, Rc and Rp as 


1 1 1 


=——+ 
Ryor _ Rs Ry Ro + Rp 


= igg(A, + A, + igg) 
Vo(A, + A) 


(22) 


Using the form of i,(t) from Eq. (20), the voltage V;(t) 
can be expressed as 


t 
ViO=Vo—Rs is + A, exp (- ) 
73 


4A, exp (- “| (23) 
1% 


Comparing the two forms of voltage V;(¢) from Eq. (19) 
and Eq. (23) gives the values of to, t; and tT, in terms of 
the known parameters as 


— CpRcRp — Ayty + 4273 


T 24 

. Ro + Rp A, + A, ad 
1 1 1)\7! 

= Goat (25) 
T 13 To 
R.C 

1) = S—M"O (Vy — isgRs) (26) 
Vor 


Now, from the expression of 73, the model parameter 
Cp can be expressed as 


Cp’ = RECué 


where 


=| 1 1 (. -)| (27) 
Tot 13 \T, 73 


The relation between Rc and Rp can be found from Eq. 
(13) and Eq. (22) as 


R C, 1 1 
D = Re( M ) (28) 
Ro + Rp tT Rs Rror — Rs 


We have shown the derivations of a two-compartment 
model based on DC measurements. Our objective is to 
express the model unknowns (Rs, Ru, Rc, Rp, Cu. Cp) 
in terms of the known parameters (iss, 73, T4, 41, A) of 
the experimentally recorded current (Eq. (20)). We are 
able to extract some parameters (Rs, Rrot, Cm; Tos T1 
T2) in this process, but the other model unknowns (Rc, 
Rp, Rm, Cp) may only be expressed in terms of some 
interdependent equations. In the next section, we will 
continue with our derivations for the parameter-extrac- 
tion of a two-compartment model, to express the ‘still- 
unknown’ model parameters in the form of independent 
equations. 


S. Pandey, M.H. White | Journal of Neuroscience Methods 120 (2002) 131-143 137 


4. Comparison with previous work 


In our derivations until now, we have avoided making 
simplifying assumptions or approximations regarding 
the relative sizes of the circuit components. At this 
point, it is interesting to review the previous work on the 
two-compartment modeling. In these works, simplifying 
assumptions are made, which provides a unique solution 
to the equivalent circuit model based on DC measure- 
ments. If we apply these assumptions to our general set 
of equations, then our equations reduce to previous 
reported work. 

Llano et al. (1991) uses the assumption of Ry = ~, 
Rp>>Rc, Rs and Cp>>Cy for modeling the 
currents of Purkinje cells. This simplifying assumption 
works well for determining the passive properties of 
Purkinje neurons. As will be shown below, the results in 
the reference (Llano et al., 1991) can be obtained from 
the general set of equations we have already derived for 
the two-compartment model. If we use the above 
assumption, then Eq. (13), Eq. (15) and Eq. (16) yield 


1 1 1 1 1 
T = CpRe, == ( + ), 
S 


T TT Cy\Ro R (29) 
T, = 73 = Cp(Ro + Rg) 
Inserting Eq. (29) into Eq. (19) provides 
Vid=V, ! Rs ( -) 
= ex 
1 fe) Ro + Rs p . 
Re t 
+ exp (30) 
Rot Rs T4 
Comparing Eq. (23) with Eq. (31) gives 
V VoR 
A, =—°— and A,=—? © (31) 
Ro + Rg RS + RsRe 


Eq. (29) and Eq. (31) are sufficient to extract the model 
parameters as 


_ 3A, _ 1%4(A, + As)” ~. Ve 
ay si AV, aa 
Oo 2%o 1 + A, (32) 
_  AVo 
. A,(A; + A) 


The results obtained by Llano et al. (1991) are the 
same as those shown in Eq. (32). Computer simulations 
are performed on a two-compartment model with the 
assumption that Rp > > Rs, Ro; Rw = 0 and Cp > > 
Cm and the results are shown in Fig. 4. A 10 mV 
command potential is applied as the voltage stimulus 
Vo. Fig. 4(A) shows the effect of series resistance 
compensation on the capacitive transient of the current 
i,(t). The model parameters chosen are Cy=5 pF, 
Cp = 100 pF, Rc = 30 MQ and Rs =5, 10 and 15 MQ. 
From Eq. (31), we see the amplitude of A> (fast) is Rc/ 


(B) 


Fig. 4. (A) The effect of series resistance compensation on the 
capacitive transient of the current i,(t) under the assumption Rp > 
> Rs, Rc; Rm = 0% and Cp> >Cm. A 10 mV command potential 
was applied as the voltage stimulus Vo. The model parameters chosen 
are Cu =5 pF, Cp = 100 pF, Rc = 30 MQ and Rg = 5, 10 and 15 MQ. 
The fitting parameters for the biexponential decay of the capacitive 
transients are (a) for Rs = 5 MQ, 13 =3.5 ms, t4 = 21.4 ps, A; = 0.28 
nA, A, =1.7 nA; (b) for Rs = 10 MQ, 73 =4 ms, t4 = 37.5 ps, A, = 
0.25 nA, A> =0.75 nA; (c) for Rs = 15 MQ, 73 =4.5 ms, t4 = 50 ps, 
A, =0.22 nA, A,=0.44 nA. (B) The effect of varying Re on the 
capacitive transient of the current i,(t) under the same assumption. 
The model parameters chosen are Cy =5 pF, Cp = 100 pF, Rs = 10 
MQ and Rc=10, 30 and 60 MQ. The fitting parameters for the 
biexponential decay of the capacitive transients are (a) for Ro =10 
MQ, 13 =2 ms, t4 = 25 ps, 4; = 0.5 nA, A. =0.5 nA; (b) for Rc = 30 
MQ, 73 =4 ms, t4 = 37.5 ps, 4; = 0.25 nA, A, =0.75 nA; (c) for Ro = 
60 MQ, 73 =7 ms, t4 = 43 ps, 4; = 0.14 nA, A> =0.86 nA. 


Rg times the amplitude of A, (slow). As the ratio of Rc/ 
Rg is increased, the amplitude of A, also increases 
(evident from Fig. 4A). The amplitude of A,;, on the 
other hand, does not change significantly. Increasing the 
series resistance compensation decreases both the fast 
(t4) and the slow (r3) time constants. Fig. 4(B) shows the 
effect of varying Rc on the capacitive transient of the 
current i,(f) under the same assumption. The model 
parameters chosen are Cy =5 pF, Cp = 100 pF, Rs = 
10 MQ and Re = 10, 30 and 60 MQ. In this case also, 
decreasing the value of Rc decreases both the fast (74) 
and the slow (13) time constants. However, the value of 
i,(0) is constant for varying values of Rc. 

Mennerick et al. (1995) use the assumption Rp > > 
Rw for modeling the passive properties of hippocampal 
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neurons. Single-exponential fits are usually inadequate 
to describe the decay of current transients from these 
hippocampal neurons, but biexponential decays pro- 
vided an adequate description of the data. The biexpo- 
nential decay of capacitive transients indicates the 
passive membrane properties may be described ade- 
quately with a two-compartment equivalent circuit 
model. As shown below, with this assumption in our 
general set of equations, we obtain the same results as 
reported in the literature (Mennerick et al., 1995). 
From Eq. (13) and Eq. (24) we have 


A,t™ + Ayt 
To = RoCp = — 4 (33 
re) c&p A; ads ) 
and with Eq. (21) and Eq. (22), 
R= Vo 1 isg(A; + Ay + iss) 
: A, + Ay + ing) Ry Vo(A, + Ay) ' (34) 
_ 1,7,(A, + A, + kel 
M Wo(Ayt4 + A573) 
Also, with Eq. (28) and Eq. (34), 
(4, + A,)(Ayt, + At)” Vo 4 
OA Ata A = ey — 
149(A; 2 4 Iss) (%4 3 (35) 
1 i 4 Ages 
Cy = : 
Ra) Apa, 


5. Two-compartment model (admittance measurements) 


A common method for measuring changes in mem- 
brane capacitances of small cells utilizes a sinusoidal 
voltage stimulus. The membrane capacitance is mea- 
sured as a function of the real and imaginary admittance 
of the cell for the case of a single-compartment model. 
In this instance, this scheme is extended for a two- 
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compartment model to help extract the remaining model 
parameters. If the AC admittance of the two-compart- 
ment model is Yror(@) = 1/Zto1r(@) = «(@)+jh(@), 
then from Fig. 3 we have 


1 
Zror(@) — Rg 
1 1 ioCyR. 
Sees ores + JOCD*p (36) 
Ry Ro + Rp + joCp Rc Rp 
which can be written as 
{a — Ro? + B’)} + JB 
{(1 — «Rs + (BRs)} 
1 R j R R 
So fp, eee iO) (37) 
Ruy R(Re + Rp) + joto) 


Equating the real parts of the Eq. (37) and using Eq. 
(13), we find 


1 
(Ro + Rp) + 22, 
a 1 1 
Ro + w13) Rg 
+ {a — Rg (0? at By} 
{(1 — as)” + (BRs)} 


Inserting Eq. (13), Eq. (27) and Eq. (38) into Eq. (28) 
gives 


Cur 


Ty 


(38) 


1 
Lat (1 +73) 


Cc 
{1S {a — Rela? + )} ] 
Rs t {(1 — ag)’ + (BRs)} 
Gates MO aga) 
To Cu Reo tT Rs Ror — Rs 


From Eq. (39), Rc can be expressed as 


ie + wto){e = Rel + Py ort tf 
1 {(l — aRs) + (BRs)} 


Re 


(40) 


{ 1 (— 1 1 ) 
ToCys \ Ty Ry — Ryor — Rs 


And from the expression C5! = R2Cyé, we have 


' 


{ 1 (“ 1 1 ) 
1 ToV Cus \ 1 Ry — Rror — Rs 


2 
Gus} 


(41) 


Cp {c + w'to){x = Re? + ay wre {Ss 
oO 


{(1 — Rs) + (BRs)"} Ty 


1 \ 1 
Rs Rror — Rs 
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Fig. 5. (A) The effect of series resistance compensation on the 
capacitive transient of the current i,(f) in a two-compartment model 
(with no approximations). The model parameters chosen are Cy = 25 
pF, Cp = 25 pF, Rc = 20 MQ, Ruy = Rp = 200 MQ and Rg = 1, 5 and 
10 MQ. The fitting parameters for the biexponential decay of the 
capacitive transients are (a) for Rs = 1 MQ, 73 =0.5 ms, t4 = 21.5 ps, 
A, =0.94 nA, 4. =9 nA, iss = 0.095 nA; (b) for Rs = 5 MQ, 73 = 0.67 
ms, T4 = 72 ps, A; = 0.7 nA, A, = 1.2 nA, igg = 0.091 nA; (c) for Rs = 
10 MQ, 73 =0.87 ms, t4= 104 ps, 4) =0.5 nA, A. =0.42 nA, iss = 
0.087 nA. (B) Plots the various values of to, tT), T2, 73 and t4 (as defined 
by Eqs. (13)—(16)) as a function of the series resistance compensation. 
The model parameters are chosen as Cy = 25 pF, Cp = 25 pF, Rc = 20 
MQ and Ru = Rp = 200 MQ. The range of the values of the other 
fitting parameters are A; =0.5—0.94 nA, A, =0.42-9.0 nA, iss = 
0.087—0.095 nA. 


Finally, we have 


1 1 1 
= : 1) and 
Ry Re \ReCusto 
1 1 1 


Ry Rror = Rs Re + Rp 


(42) 


From the results obtained in previous sections, all the 
six two-compartment model parameters (Rc, Rp, Rm, 
Rs, Cy, Cp) are independently expressed in terms of the 
known parameters (iss, 13, 74, 41, 42) of the current i;(t) 
expression (Eq. (20)) and the known parameters («(@), 
B(@), @) of the admittance measurements. A flowchart 
summarizing the steps for the extraction of the model 
parameters is shown at the end of the next section as 


nA 


mS = time 
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Fig. 6. (A) The effect of varying Rc on the capacitive transient of the 
current i,(¢) in a two-compartment model (with no approximations). 
The model parameters chosen are Cy =25 pF, Cp =25 pF, Rs =5 
MQ, Ry = Rp = 200 MQ and Rc =2, 20 and 100 MQ. The fitting 
parameters for the biexponential decay of the capacitive transients are 
(a) for Ro =2 MQ, 13 =0.3 ms, t4 = 19.3 ps, A) = 1.7 nA, Ar = 0.22 
nA, iss = 0.095 nA; (b) for Rc = 20 MQ, 73 = 0.67 ms, t4 = 72 ps, Ay = 
0.69 nA, A> = 1.22 nA, isg = 0.09 nA; (c) for Rc = 100 MQ, 73 = 1.84 
ms, t4 = 103 ps, Ay = 0.19 nA, A> = 1.73 nA, igs = 0.08 nA. (B) Plots 
the various values of to, 1), T2, T3 and ty (as defined by Eqs. (13)—(16)) 
as a function of Rc. The model parameters are chosen as Cy = 25 pF, 
Cp = 25 pF, Rs =5 MQ and Ry = Rp = 200 MQ. The range of the 
values of the other fitting parameters 4; = 0.12—1.7 nA, Az = 0.22-1.8 
nA, iss = 0.073—0.095 nA. 


Fig. 11. Such an algorithm can be easily implemented in 
software, which takes (as input) the DC and AC 
measurement results and gives (as output) the values 
of the model parameters. 


6. Two-compartment model (simulation results) 


Computer simulations of the two-compartment 
model, based on our general parameter-extraction 
method, were performed on XSPICE. The effect of 
varying the series resistance compensation on the two- 
compartment model current response is shown in Fig. 5. 
The model parameters chosen are Cy = 25 pF, Cp = 25 
pF, Rc = 20 MQ, Ruy = Rp = 200 MQ and Rg = 1, 5 and 
10 MQ. Fig. 5(A) shows the effect of series resistance 
compensation on the capacitive transient of the current 
i,(t) in a two-compartment model. From Rs = 1 MOQ to 
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Fig. 7. (A) The effect of varying Rp on the capacitive transient of the 
current i,(¢) in a two-compartment model (with no approximations). 
The model parameters chosen are Cy =25 pF, Cp =25 pF, Rs =5 
MQ, Ru = 200 MQ, Rc = 20 MQ and Rp = 20, 100 and 400 MQ. The 
fitting parameters for the biexponential decay of the capacitive 
transients are (a) for Rp =20 MQ, 13=0.39 ms, t4=60 pts, A; = 
0.725 nA, A. =1.02 nA, iss = 0.26 nA; (b) for Rp = 100 MQ, 73 = 
0.615 ms, t4 = 70.3 ps, A; = 0.68 nA, Az = 1.2 nA, igg = 0.125 nA; (c) 
for Rp = 400 MQ, 73 = 0.71 ms, t4 = 73 ps, Ay = 0.7 nA, Az = 1.23 nA, 
igs = 0.07 nA. (B) Plots the various values of to, tT), T2, 73 and T4 (as 
defined by Eqs. (13)—(16)) as a function of Rp. The model parameters 
are chosen as Cy =25 pF, Cp =25 pF, Rs =5 MQ and Ry = 200 
MQ. The range of the values of the other fitting parameters 4; = 0.7— 
0.76 nA, Az = 0.95—1.23 nA, igg = 0.07—0.3 nA. 


Rs=10 MQ, the increase in 13 (slow) is not so 
significant (0.5—0.87 ms) compared to the increase in 
t4 (fast) (21.5—104 ps). Also, the value of i)(0) is very 
sensitive to changes in Rs. Fig. 5(B) plots the various 
values of to, 7, T2, 73 and ty (as defined by Eqs. (13)- 
(16)) as a function of the series resistance compensation. 
According to Eq. (13), to is relatively independent of 
Rs, as shown in Fig. 5. It should be noted here that fast 
exponential decay time constant of the capacitive 
transient is t4, while the slow exponential decay time 
constant is 13. The variables t, (fast) and t (slow) 
represent the two time constants for a two-compartment 
model under the approximation of Ry =o, Rp > > 
Rc, Rs and Cp > > Cy (Llano et al., 1991). As seen 
in Fig. 5(B), the time constants t,; and ty are close to 
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Fig. 8. (A) The effect of varying Ry on the capacitive transient of the 
current i,(¢) in a two-compartment model (with no approximations). 
The model parameters chosen are Cy =25 pF, Cp =25 pF, Rs =5 
MQ, Rp = 200 MQ, Re = 20 MQ and Ry = 10, 50 and 200 MQ. The 
fitting parameters for the biexponential decay of the capacitive 
transients are (a) for Ry = 10 MQ, 13 =0.6 ms, t4 = 56 ps, A; = 0.46 
nA, A,=1.05 nA, iss = 0.7 nA; (b) for Rw =50 MQ, 73 =0.66 ms, 
T4 = 69 us, A; = 0.62 nA, A, = 1.2 nA, igg = 0.22 nA; (c) for Ry = 200 
MQ, 13 = 0.67 ms, t4=72 ws, 4; =0.7 nA, A. =1.22 nA, iss = 0.09 
nA. (B) Plots the various values of to, |, 12, 73 and ty (as defined by 
Eqs. (13)—(16)) as a function of Ry. The model parameters are chosen 
as Cy = 25 pF, Cp = 25 pF, Rs = 5 MQ and Rp = 200 MQ. The range 
of the values of the other fitting parameters A; = 0.44—0.71 nA, 42 = 
1.03-1.23 nA, iss = 0.07—0.53 nA. 


each other, but the difference between the time constants 
T and 713 increases with increasing values of Rs. This 
shows the limitation of using an approximated two- 
compartment model with a high series resistance. In 
addition, the difference between the time constants 1 
and t3 increases with an increasing value of the 
connecting resistance Rc. Finally, both time constants 
of the general two-compartment model, 73 and 14, 
depend and change with varying values of Rg. This is 
in contrast to the finite-cable model, where changing Rs 
affects only the first, faster time constant and leaves the 
slower one unchanged. 

Fig. 6 shows the effect of varying Rc on the current 
response of a two-compartment model. Fig. 6(A) shows 
the effect of varying Rc on the capacitive transient of the 
current i;(¢). The model parameters chosen are Cy = 25 
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Fig. 9. (A) The effect of varying Cy, on the capacitive transient of the 
current i,(¢) in a two-compartment model (with no approximations). 
The model parameters chosen are Cp =25 pF, Rs =5 MQ, Ry = 
Rp = 200 MQ, Rc = 20 MQ and Cy =S, 25 and 100 pF. The fitting 
parameters for the biexponential decay of the capacitive transients are 
(a) for Cu = 5 pF, 73 = 0.6 ms, t4 = 18 pts, 4) = 0.42 nA, A. = 1.5 nA, 
igg = 0.09 nA; (b) for Cy =25 pF, 13 = 0.67 ms, t4 = 72 ps, A; =0.7 
nA, A> = 1.22 nA, iss = 0.09 nA; (c) for Cy = 100 pF, 73 = 1.03 ms, 
T4 = 175 ps, Ay = 1.3 nA, A, =0.63 nA, iss = 0.09 nA. (B) Plots the 
various values of To, 71, T2, 73 and 74 (as defined by Eqs. (13)—(16)) asa 
function of Cy. The model parameters are chosen as Cp =25 pF, 
Rs =5 MQ and Ru = Rp = 200 MQ. The range of the values of the 
other fitting parameters are A, =0.42—-1.6 nA, A, =0.35-1.5 nA, 
iss = 0.09 nA. 


pF, Cp =25 pF, Rs =5 MQ, Ry = Rp = 200 MQ and 
Rco=2, 20 and 100 MQ. Fig. 6(B) plots the various 
values of to, T, T2, 73 and ty (as defined by Eqs. (13)- 
(16)) as a function of Re. The model parameters chosen 
in this case are the same as those for Fig. 6(A). From 
both the plots, the slow time constant 73 is very sensitive 
to increases in the value of Rc, whereas the other time 
constant ty is relatively constant with any changes in the 
value of Rc. Also, the difference between 1. and 73 
becomes smaller with increasing values of Rc. This 
shows that the approximation of Cy =0 by Nadeau and 
Lester (2000) is valid only for larger values of Rc where 
the difference between the two time constants, 73 and t4, 
is large. 

Fig. 7 shows the effect of varying Rp on the current 
response of a two-compartment model. The model 
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Fig. 10. (A) The effect of varying Cp on the capacitive transient of the 
current i,(¢) in a two-compartment model (with no approximations). 
The model parameters chosen are Cy =25 pF, Rs =5 MQ, Ru = 
Rp = 200 MQ, Re = 20 MQ and Cp =5, 25 and 100 pF. The fitting 
parameters for the biexponential decay of the capacitive transients are 
(a) for Cp = 5 pF, 13 = 0.23 ms, t4 = 39 ps, A; = 1.4 nA, A, =0.52 nA, 
iss = 0.09 nA; (b) for Cp = 25 pF, 13 =0.67 ms, t4 = 72 ps, A; =0.7 
nA, A, = 1.22 nA, iss = 0.09 nA; (c) for Cp = 100 pF, 73 = 2.33 ms, 
T4 = 90 ps, A; = 0.44 nA, A. = 1.47 nA, iss =0.09 nA. (B) Plots the 
various values of to, 7), T2, 73 and ty (as defined by Eqs. (13)—(16)) as a 
function of Cp. The model parameters are chosen as Cy = 25 pF, 
Rs =5 MQ, Ry = Rp = 200 MQ and Rc =20 MQ. The range of the 
values of the other fitting parameters are 4; =0.39-1.39 nA, 4. = 
0.52—1.52 nA, igg = 0.09 nA. 


parameters chosen are Cy = 25 pF, Cp =25 pF, Rs =5 
MQ, Ru = 200 MQ, Rc = 20 MQ and Rp = 20, 100 and 
400 MQ. Fig. 7(A) shows the effect of varying Rp on the 
capacitive transient of the current 7\(¢). Fig. 7(B) plots 
the various values of to, T), T2, T3 and ty as a function of 
Rp, choosing the same values of the model parameters. 
Similar to Fig. 6, we see a strong dependence of the slow 
time constant t3 on the varying values of Rp, while the 
fast time constant t4 is relatively independent of such 
changes in Rp. The difference between 12 and 73 is more 
pronounced for smaller values of Rc. 

Fig. 8 shows the effect of varying Ry, on the current 
response of a two-compartment model. The model 
parameters chosen are Cy =25 pF, Cp =25 pF, Rs = 
5 MQ, Rp = 200 MQ, Rc = 20 MQ and Ry = 10, 50 and 
200 MQ. Fig. 8(A) shows how the capacitive transient of 


142 S. Pandey, M.H. White | Journal of Neuroscience Methods 120 (2002) 131-143 


Perform DC Measurements, 
Record the current i,(t) 


Use Curve-Fitting Method to 
find the parameters of eq.(20) 
Liss, T3, Ta, Ai, Aa] 


Determine the values of Rror, 
Rg, Cy using eq.(21) 


Determine the values of these 
unknowns: To, 7), T2, & 


Perform AC Measurements 
Find a(@) and B(@) 


Determine values of Rc, Cp, 
Rp, Ru using eqs.(40-42) 


Fig. 11. A flowchart to summarize the steps for the parameter- 
extraction of a two-compartment model using the algorithm described 
in this work. 


the current 7,(¢) varies with changing Ry. Fig. 8(B) plots 
the various values of to, T), T2, T3 and ty as a function of 
Ry, choosing the same values of the model parameters. 
As seen in Fig. 8, the time constants t,, t2, 73 and t4 
exhibit less dependence on the varying values of Ry. 
Yet, there is significant difference between the curves of 
T2 and 13, suggesting that for relatively smaller values of 
Rc, the more accurate time constant (slow) is 73 and not 
T2 (which is the slow time constant under the approx- 
imation made by Llano et al. (1991)). 

Fig. 9 shows the effect of varying Cy, on the current 
response of a two-compartment model. The model 
parameters chosen are Cp = 25 pF, Rg =5 MQ, Ry = 
Rp = 200 MQ, Rc = 20 MQ and Cy =5, 25 and 100 pF. 
Fig. 9(A) shows the effect of varying Cy on the 
capacitive transient of the current i,(¢). Fig. 9(B) plots 
the various values of to, 7), 72, 73 and t4 as a function of 
Cm, choosing the same values of the model parameters. 
As seen in Fig. 9, both the time constants 73 and t4 vary 
with changing values of Cy; the dependence of the slow 
time constant 13 on such changes being much more 
significant. The curve of 1, is relatively independent of 
varying Cy because it was the slow time constant for the 


case where Cy, = 0 in a two-compartment model (Llano 
et al., 1991; Nadeau and Lester, 2000). In this case, there 
is noticeable error involved in the computation of the 
time constants if approximations are made in the two- 
compartment model. This error becomes more and more 
significant as the value of Cy, increases. 

Fig. 10 shows the effect of varying Cp on the current 
response of a two-compartment model. The model 
parameters chosen are Cy = 25 pF, Rs =5 MQ, Ruy = 
Rp = 200 MQ, Rc = 20 MQ and Cp = 5, 25 and 100 pF. 
Fig. 10(A) shows the effect of varying Cp on the 
capacitive transient of the current i)(¢). Fig. 10(B) plots 
the various values of to, tT), T2, T3 and t4 as a function of 
Cp, choosing the same values of the model parameters. 
As seen from Fig. 10, the fast time constant t4 is 
relatively constant with changing values of Cp, while the 
slow time constant 13 shows a sharp, linear increase with 
increasing values of Cp. The curves of 12 and 13 follow 
each other closely and so does the curves of t; and t4. As 
such, making an approximation about the relative size 
of Cp does not introduce any significant error in the 
computation of the time constants. 


7. Conclusion 


In summary, we have derived a general solution for 
the parameter-estimation of a two-compartment model. 
With simplifying approximations, used by previous 
authors, our derivations are consistent with their results. 
The computer simulation results of the general two- 
compartment model give valuable insight into the role of 
each model parameter in the current response of such an 
equivalent circuit. We have compared the various time 
constants (those from approximated models and those 
from our general model) under the effects of varying 
model parameters. The difference between the time 
constants (related to the approximated models and our 
general model) increases with decreasing values of the 
connecting resistance Rc. As seen from the simulations, 
each of the model parameters has a unique effect on the 
variation of the different time constants. Employing an 
approximated model may not reveal the exact values of 
the model parameters, which are calculated from these 
time constants. The general method for the extraction of 
parameters is not limited by any simplifying approxima- 
tions and so gives an exact solution of the two- 
compartment model. 
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